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' We present a statistical mechanical model of aggregation in colloidal systems with DNA mediated 

I interactions. We obtain a general result for the two-particle binding energy in terms of the 

, hybridization free energy AG of DNA and two model dependent properties: the average number 

of available DNA bridges (A'^) and the effective DNA conccentration Ce//. We calculate these 
parameters for a particular DNA bridging scheme. The fraction of all the n-mers, including the 
infinite aggregate, are shown to be universal functions of a single parameter directly related to the 
two-particle binding energy. We explicitly take into account the partial ergodicity of the problem 
' resulting from the slow DNA binding-unbinding dynamics, and introduce the concept of angular 

04 , localization of DNA linkers. In this way, we obtain a direct link between DNA thermodynamics 

and the global aggregation and melting properties in DNA-coUoidal systems. The results of the 
theory are shown to be in quantitative agreement with two recent experiments with particles of 
micron and nanometer size. PACS numbers: 81.16.Dn, 82.20.Db, 68.65.-k, 87.14. Gg 
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I. INTRODUCTION 



In the past ten years, there have been a number of advances in experimental assembly of nanoparticles with 
DNA-mediated interactions ([1], [2], [3], [4], [5], [6]). While this approach has a potential of generating highly organized 
^ and sophisticated structures([7],[8]), most of the studies report random aggregation of colloidal particles([9],[10]). 
I , Despite these shortcomings, the aggregation and melting properties may provide important information for future 
' O ■ development of DNA-based self-assembly techniques. These results also have more immediate implications. For 
instance, the observed sharp melting transition is of particular interest for biosensor applications [11]. For these 
reasons the development of a statistical mechanical description of these types of systems is of great importance. It 
should be noted that the previous models of aggregation in coUoidal-DNA systems were either phenomenological or 
oversimplified lattice models([12],[13],[14]), which gave only limited insight into the physics of the phenomena. 
' In this paper, we develop a theory of reversible aggregation and melting in coUoidal-DNA systems, starting from 
, the known thermodynamic parameters of DNA (i.e. hybridization free energy AG), and geometric properties of 
DNA-particle complexes. The output of our theory is the relative abundance of the various colloidal structures 
, formed (dimers, trimers, etc.) as a function of temperature, as well as the temperature at which a transition to an 
\0 ' infinite aggregate occurs. The theory provides a direct link between DNA microscopies and experimentally observed 
' morphological and thermal properties of the system. It should be noted that the hybridization free energy AG 
^5 depends not only on the DNA nucleotide sequence, but also on the salt concentration and the concentration of DNA 
linker strands tethered on the particle surface[15]. In this paper AG values refer to hybridization between DNA free 
in solution. 

In a generic experimental setup, particles are grafted with DNA linker sequences which determine the particle 
^ [ type(^ or B). In this paper we will restrict our attention to a binary system^. These linkers may be flexible or 
■ rigid. A selective, attractive potential between particles of type A and B can then be turned on by joining the linkers 
C5 [ to form a DNA bridge. This DNA bridge can be constructed directly if the particle linker sequences are chosen to 
O have complementary ends. Alternatively, the DNA bridge can be constructed with the help of an additional linker 
^ ' DNA. This additional linker is designed to have one end sequence complementary to the linker sequence of type A 
J> . particles, and the other end complementary to type B. The properties of the DNA bridge formed will depend on the 

hybridization scheme (see figure 1). 
^ , The plan for the paper is as follows. In section II. A we provide a description of the problem. In section II. B 
' we determine the bridging probability for the formation of a DNA bridge between two colloids, assuming the known 
thermodynamic parameters of DNA(hybridization free energy AG). Using this bridging probability as input, in 
section II. C we calculate the effective binding free energy e^s for the formation of a dimer. Sections III.A-III.C 
establish the connection between the theory and the experimentally determined melting profile /(T), the fraction of 



This restriction to binary systems is consistent with the current experimental approach. In a recent work we demonstrated that if each 
particle has a unique linker sequence, one might be able to programmably self-assemble nanoparticle clusters of desired geometry [8]. 
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FIG. 1: (Color Online) Graphical depiction of various schemes for DNA bridging. A) A freely-jointed, rigid bridge constructed 
from complementary linker DNA. B) A flexible bridge can be constructed using complementary linker DNA. C) A rigid bridge 
constructed from short, flexible linker DNA and a long, rigid linker. 

unbound particles as a function of temperature. In particular, we demonstrate how knowledge of eAB can be used to 
determine this profile, including the effects of particle aggregation. In Section III.D the theory is compared with two 
recent experiments detailing the reversible aggregation of colloids with DNA-mediated attraction [12], [9]. The main 
results of the model are summarized in section IV. 



II. DNA-MEDIATED INTERACTIONS 

A. Description of the problem 

We consider particles of type A and B which form reversible AB bonds as a result of DNA hybridization. The 
task at hand is to determine the relative abundance of the various colloidal structures that form as a function of 
temperature. From this information we can determine which factors affect the melting and aggregation properties in 
DNA-colloidal assemblies. To do so we must determine the binding free energy for all of the possible phases(monomer, 
dimer, infinite aggregate), and then apply the rules for thermodynamic equilibrium. As we will see in section III, 
these binding free energies can all be simply related to eAB, the binding free energy for the formation of a dimer. Our 
task is thus reduced to determining ^ab from the thermodynamic parameters of DNA and structural properties of the 
DNA linkers. In our statistical mechanical framework, eAB is calculated from the model partition function, taking 
into account the appropriate ensemble averaging for the non-ergodic degrees of freedom. The result is related to the 
bridging probability for a pair of linkers. By considering the specific properties of the DNA bridge that forms, the 
bridging probability can be related to the hybridization free energy AG of the DNA. In this way, we obtain a direct 
link between DNA thermodynamics and the global aggregation and melting properties in coUoidal-DNA systems. 

B. Bridging probability 

To begin we relate the hybridization free energy AG for the DNA in solution to the bridging probability for a pair 
of linkers. This bridging probability is defined as the ratio ''^J""'' , with Pbound the probability that the pair of linkers 
have hybridized to form a DNA bridge, and Pfree the probability that they are unbound. This ratio is directly 
related to the free energy difference of the bound and unbound states of the linkers AG (throughout this paper we 
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will use units with ks = 1): 



bound 



exp 



free 



Ceff 



-AG 



Ceff 



exp 



-AG 



J P(ri,r)P(r2,r)d3r 
(/P(r,r')d3r)' 



(1) 
(2) 



Here Co = IM is a reference concentration. P(r, r') is the probability distribution function for the linker chain which 
starts at r' and ends at r. The effective concentration Ce// is a measure of the change in conformational entropy 
of the linker DNA as a result of hybridization. It will depend on the properties of the linker DNA(ex: flexible vs. 
rigid), and the scheme for DNA bridging(ex: hybridization of complementary ends vs. hybridization mediated by an 
additional linker). Ce// is the concentration of free DNA which would have the same hybridization probability as the 
grafted linkers in our problem. As discussed in section III.D. the DNA linker grafting density also plays an important 
role in determining the possible linker configurations and hence Ce/ / . 

Assuming that the size of the linkers is much smaller than the particle radius R, we first consider the problem 
in a planar geometry. Let the two linkers be attached to two parallel planar surfaces separated by a distance 2h. 
Referring to figure 2 we see that r' is the location where the linker DNA is grafted onto the particle surface, and r is 
the position of the free end. 




FIG. 2: (Color Online) The statistical weight of a bound state is calculated by determining the number of hybridized configu- 
rations for two complementary linker chains relative to the number of unhybridized configurations. 

In this paper we consider hybridization by complementary, rigid linker DNA (scheme A in Figure 1). This scheme 
is particularly interesting since it is directly related to several recent experiments [9], [12]. In a future work we will 
address other hybridization schemes. We assume that L < Lp and L <^ R, where Lp ~ 50nm is the persistence 
length of ds DNA and L is the ds linker DNA length. In this regime, the linker chains can be treated as rigid rods 
tethered on a planar surface. The interaction is assumed to be point-like, in which a small fraction A/L of the linker 
bases hybridize. 

We can calculate the effective concentration by noting that the overlap integral in Eq. (2) is proportional to the 
volume of intersection of two spherical shells (red and blue circles in Figure 3) : 



Ce// 



2TTrA e (L - |r; - /2) 



{27tL^AY 



27r2L2 



(3) 



here A = A^/sin/? and r = y _ _ i-'2| /4 (see notations in Figure 3). We have used the fact that cos/3/2 = 
k'l ~ \ /2i. Ceff and the binding probability are largest when the linkers are grafted right in front of each other. 
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i.e. when |r'j^ — ~ By taking the limit /i « L we arrive at the following result for the corresponding "bridging" 
free energy 

AG A ~ AG A + T log [4ttL^Co] . (4) 

This free energy remains nearly constant for any pair of linkers, as long as they can be connected in principle, i.e. 
\r'i — r'jl < 2L. This limits the maximum lateral displacement of the linkers: r± < — /i^, and therefore sets the 

effective cross-section of the interaction: 

a= Trri = 47r (L^ - h^) (5) 




FIG. 3: (Color Online) Cross-sectional view of the hybridization of two complementary rigid linker DNA. The effective 
concentration is calculated in a planar approximation to the particle surface. 



C. Effective binding free energy 



We now proceed with the calculation of the effective free energy eAB , which is associated with the formation 
of a dimer from a pair of free particles, A and B. Since the DNA coverage on the particle surface is not uniform, 
this free energy, and the corresponding partition function Z, would in principle depend on the orientations of the 
particles with respect to the line connecting their centers. The equilibrium binding free energy would correspond to 
the canonical ensemble of all possible orientations, i.e. sab = — T'log47r {Z). However, this equilibrium can only be 
achieved after a very long time, when the particle pair samples all possible binding configurations, or at least their 
representative subset. The real situation is different. After the first DNA-mediated bridge is created the particle 
pair can still explore the configurational space by rotating about this contact point. However, after the formation of 
two or more DNA bridges (at certain relative orientation of the particles), the further exploration requires multiple 
breaking and reconnecting of the DNA links, which is a very slow process. We conclude that the system is ergodic 
with respect to the various conformations of the linker DNA for fixed orientations of the particles, but the orientations 
themselves are non-ergodic variables. The only exceptions are the single-bridge states: the system quickly relaxes to 
a more favorable orientational state (unless the DNA coverage is extremely low, and finding a second contact is very 
hard). If N denotes the number of DNA bridges constituting the AB bond, the appropriate expression for eAB hi 
this partially ergodic regime is the so-called component averaged free energy([16],[17]): 



<^AB = -T(logZ)^>2 



(6) 
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Each DNA bridge between particles can be either open or closed. 



^bridge — 1 + exp 



(7) 



Here r- is the 2D position where the bridge is grafted onto surface i. We now consider a generic case when the 
interaction free energy AG depends on the separation between planar surfaces 2h' , and the separation of grafting 
points v'l — r'2, without assumption of a particular bridging scheme. If the probability for bridge formation is small, 
two DNA linkers on the same surface will not compete for complementary linkers. In this regime the free energy can 
be calculated by summing over the contribution from each bridge that forms between dimers. 



1 + exp 



AG(/;'.r:-r') 



(8) 



We convert the summation to integration by introducing the linker areal grafting density a 

AG(/i',r'i-r^) 



^ = -^1 J ^i(r'i)a2(r^)log 



1 + exp 



(9) 



Changing variables to Ar = — and p = (r'^ + rg) /2, we can reintroduce the notion of a bridging cross-section 
a{h'), this time in a model-independent manner: 



a{h')log 



1 + exp 



AGo{h') 



I 



d^Arlog 



1 + exp 



AG(/i',Ar) 



(10) 



Here AGo{h') = AG{h', Ar = 0) is the minimum free energy with respect to the separation between grafting points 
Ar. We can now write the free energy: 



P = -^ j (ri{-p)a2{'p)a{h')log 



1 + exp 



AGoih') 



(11) 



We now convert from the planar geometry to the spherical particle geometry using the Derjaguin approximation [18]. 



d p = pdpdcj) 



h' = h 



2R 



(12) 
(13) 



Let AG* be the minimal value of the bridging free energy. Then the result for F can be rewritten as: 



F = -TN log 



1 + exp 



AG* 



(14) 



Here A'' has a physical meaning as the number of potential bridges for given relative positions and orientations of the 
particles: 



N = J aiC^)a2(p)a(/i') 



log 



l + exp(-AGo(/i')/T 



log 



1 + exp 



(-AG*/r) 



d^p 



(15) 



One can calculate the average value of N in terms of the average grafting density, a = (cti) = (cr2) : 



(AT) = 27rRa' 



■Ja{h') 



log 


1 + exp (-AGo(/i')/T 


)] 


log 


1 + exp (-AGJT 







dh' 



(16) 
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In a generic case of randomly grafted linkers, (N) completely defines the overall distribution function of N, which 

must have a Poisson form: P{N) = - — '-j^, . The average number of bridges (TV) between two particles depends 

on both the DNA linker grafting density a and the bridging probability determined from AG. 

The free energy for the formation of a dimer eab = {P)2+ ~ Tlogil. The second term is the entropic contribution 
to the free energy, which comes from integration over the orientational and translational degrees of freedom of the 
second particle. Because the system is not ergodic in these degrees of freedom, the accessible phase space ft will be 
reduced by a factor of P2-1-. P2+ is the probability that there are at least two DNA bridges between the particles. 
In terms of the average number of bridges (A^") between particles, we have the following relations: 

P2+ = l-(l + (17) 

^ <^^K^ ^^^^ 



eAB = -T^(Ar)2+log 



1 + exp — - 



log [P2+4nS{2Rfco\ } (19) 



Here S is the localization length of the AB bond, which comes from integrating the partition function over the radial 
distance between particles. 

We now can calculate {N) for the case of freely-jointed rigid bridging considered earlier (i.e. for scheme A). In a 
previous section we provided a direct calculation of the interaction free energy, AGo (h) « const = AG* (eq.4), and 
bridging cross-section, a{h') = 47r(L^ — h'^). Applying eq. 16 we arrive immediately at the following result. 

L 

{N) = Sn^a^R J [L^ - h'^)dh' = l^^'^'^^' (20) 


Having determined the free energy, we are now in a position to determine the melting properties for DNA colloidal 
assemblies. 



III. AGGREGATION AND MELTING IN COLLOIDAL-DNA SYSTEMS 



At this stage of the paper we have calculated the binding free energy eab for an AB pair, starting with the 
thermodynamic parameters of DNA (hybridization free energy AG). In this section of the paper we establish the 
connection between that result and the experimentally observable morphological behavior of a large system. One of 
the ways to characterize the system is to study its melting profile /(T), which is the fraction of imboimd particles 
as a function of temperature. To determine the profile we calculate the chemical potential for each phase(monomer, 
dimer, etc.) and apply the thermodynamic rules for phase equilibrium. We will demonstrate how the single binding 
free energy eAB can be used to determine the contribution of each phase to the melting profile, including the effects 
of aggregation. 



A. Dimer formation 



To begin we discuss the formation of dimers via the reaction A + B ^ 
of the i^^ species in terms of the particle concentrations Ci = 

II A = Tlog (ca) 
Mb =2" log (cs) 
Mas = ^log {cab) + eAB 



AB. We can express the chemical potential 

(21) 
(22) 
(23) 



Here cab is the binding free energy for the formation of a dimer. In terms of the potential V{r) between A and B 
type particles we have: 



eAB 



-Tlog 



4n{2Rfco j drexp 



V{r) 



(24) 
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In this section we are not particularly concerned with the specific form of the DNA- induced potential V{r), having 
already determined ej^g in the previous section. We simply note that the prefactor 47r(2i?)^ arises since the interaction 
is assumed to be isotropic, with R the particle radius. Equilibrating the chemical potential of the various particle 
species, we obtain the condition for chemical equilibrium. 



Ma + Mb — Mab 

The result is a relationship between the concentration of dimers and monomers. 

caCb \-€ab 
Cab = exp 



(25) 



(26) 



The overall concentration of particles in monomers and dimers must not differ from the initial concentration. 

C\ =CA + CAB (27) 
4 = CB + CAB (28) 



'B 



If the system is prepared at equal concentration, c\ 
CA = cb = c. Written in terms of the fraction of unbound particles / 
unbound fraction. 



ictot, subtracting the two equations we see that 
- — we have a quadratic equation for the 



1 = / + exp 



-CAB 



To simplify we have defined an effective free energy cab for the formation of a dimer. 



CAB = CAB - T log 



Ctot 



2Cn 



(29) 



(30) 



The solution for the fraction of unbound particles as a function of temperature is simply: 









-1 + Y'l + 4exp 


— «AB 
T 


2 exp 


— EAB 
T 





(31) 



Previous studies [9] only included the dimer contribution to the melting properties of DNA colloidal assemblies. 
With the basic formalism at hand, we can now extend the preceding analysis to include the contribution of trimers 
and tetramers. 



B. Trimers and tetramers 



Now consider the formation of a trimer via 2 A + B ^ ABA. The chemical potential is slightly different in this 
case. 

Maba = ^log {cab a) + eABA (32) 

Taking into account that there are now two AB bonds in the structure, one might conclude that cab a = 2eAB- This 

is not quite correct, since there is a reduction in solid angle available to the third particle. To form a trimer, an AB 
bond forms first, which contributes cab to cab A- Some simple geometry shows that the remaining A particle only 
has 37r steradians of possible bonding sites to particle B. Making this change in the prefactor of eq. 24, one can see 
that the second bond contributes cab — T\og (|) to caba- 

e AB A = 2eAB - T log (^^^ (33) 

The equation for chemical equilibrium can once again be expressed in terms of the particle concentrations. 

2Ma + Mb = Maba (34) 

3c^CB r-2eABl 
CABA = -^—^exp — - — (35) 
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To include the trimer contribution, we note that there are two possible varieties, with cab a = ^bab- 



c\ = CA + CAB + 2CAByl + CbAB 
Cb = Cb + CAB + CABA + '^CbAB 

Following the same line of reasoning as before, the resulting equation for the unbound fraction / is: 



1 = / + exp 



-^AB 



-exp 



-2? 



AB 



(36) 
(37) 



(38) 



For tetramers we will follow the same general reasoning, however in this case there are two different structure types. 
The reaction 2A + 2B ^ ABAB results in the formation of string like structures. 



IJ-ABAB = Tlog {caBAb) + CABAB 



(39) 



As in the trimer case, the last particle has Stt steradians of possible bonding sites, and contributes cab — Tlog (|) to 

^ABAB- 



CABAB = 3eAB - T log 

+ = t^ABAB 



CABAB 



Q\ 2 2 2 
O \ Cj^Cg 



exp 



-Scab 



(40) 

(41) 
(42) 



If an A type particle approaches a trimer of variety ABA, a branched structure can result. The reaction ^A + B 
AAAB results in the formation of these branched structures. 



I^AAAB = Tlog (caAAb) + ^AAAB 



(43) 



For the branched case, the last particle has approximately 2'7r steradians of possible bonding sites, and contributes 

eAB - Tlog (i) to CAAAB- 



iAAAB = Scab - T log - 



3/iA + Mb — ^^AAAB 



caaab 



CaCB 



exp 



-3eAi 



(44) 
(45) 
(46) 



To include all of the tetramer contributions, note that there are two branched varieties, with caaab = cbbba- Finally 
we impose the constraint that the initial particle concentrations do not differ from the concentration of all the n-mers, 
for n=l,2,3,4. 



c'a = ca + Cab + "^caba + cbab + "^cabab + ^caaab + cbbba 
Cb = Cb + Cab + cab a + "^cbab + '^cabab + caaab + ^cbbba 



(47) 
(48) 



The final result is an equation for the unbound fraction / expressed entirely in terms of the effective free energy cab 
of a dimer. 



1 = / + exp 



-CAB 



T 



/' + Texp 



-2?Ai 



T 



21 



■ exp 



-re. 



AB 



T 



(49) 



For high temperatures, the melting profile is governed by the solution to this polynomial equation for /. For 
temperatures below the melting point we expect to find particles in large extended clusters. We now proceed to 
calculate the equilibrium condition between monomers in solution and the aggregate. 
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(-ab/T 



FIG. 4: (Color Online) The actual unbound fraction / is the concatenation of the aggregate profile for T <T* and the n-mer 
profile for T > T* . The fraction of particles in dimers, trimers, and tetramers is also plotted. 



C. Reversible sol-gel transition 

To understand the basic structure of the aggregate, we simply note that there are many DNA attached to each 
particle. This gives rise to branching, as in the discussion of possible tetramer structures. Since the DNA which 
mediate the interaction are grafted onto the particle surface, once two particles are bound, the relative orientation of 
the pair is essentially fixed. The resulting aggregate is a tree-like structure, and the transition to an infinite aggregate 
at low temperatures is analogous to the sol-gel transition in branched polymers[19]. 

Particles in the aggregate are pinned down by their nearest neighbor bonds, so we do not consider their translational 
entropy. As a result the chemical potential is simply = too- Equilibrating the chemical potential of the monomer 
in solution and in the aggregate we have: 

riog(c)=eoo (50) 
Eoo = e^B - riog(7^) (51) 

Here ~ 1 is the configurational entropy of the branched aggregate, per particle. 

The concentration of particles in the aggregate Cqo is the the total concentration minus the n-mer concentration. 
Here ci = c^i -I- is the total monomer concentration, C2 — cab is the total dimer concentration, etc. 

Coo ~ Ctot - Ci - C2 - C3 - C4 (52) 

Expressed in terms oilAB and the fraction of solid angle available to particles in the aggregate — have: 

CAB 



, 1 
/oo = exp 



T 



(53) 



The transition from dimers, trimers, etc. to the aggregation behavior is the temperature T* at which foo{T*) is a 
solution to eq. 49. In words, T* is the temperature at which the aggregate has a non-zero volume fraction. The 
fraction of unbound particles for these colloidal assemblies will be governed by eq. 53 for T < T* and eq. 49 for 
T > T* . As claimed, we can simply relate the unbound fraction toT ab for both n-mers and the aggregate. 



D. Comparison to the experiment 



Let's consider the experimental scheme of Chaikin et al[9]. In the experiment, R = .5fim polystyrene beads were 
grafted with ds DNA linkers of length L ~ 20nm. The 11 end bases of the A and B type particles were single stranded 
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and complementary. We have already determined the bridging probability in this scenario(see section TLB). In the 
experiment [9] a polymer brush is also grafted onto the particle surface, which will have the effect of preferentially 
orienting the rods normal to the surface(See Figure 3). This confinement of the linker DNA can be incorporated 
quite easily into our results for AG and (N). Returning to section II. B, when integrating over linker conformations 
we simply confine each rigid rod to a cone of opening angle 2a. The upper bound for the polar integration is now a 
as opposed to tt. 



AG A ^ AG A + T\og [4:TtL^Co{1 - cosaf] 



(54) 



The alignment effect should also be taken into account when calculating (N) . If the particles are separated by less 
than 2Lcosa the end sequences will be unable to hybridize. Following the same steps as before, the lower bound for 
the h integration is now Lcosa as opposed to 0. 



{L^ - h'^)dh 



L cos a 
3 



1 



COS a 



(cos^ a — 3) 



(55) 
(56) 



In the absence of the brush, and at sufficiently low linker grafting density cr, the alignment effect could be removed 
by setting a = ^, in which case we recover our previous results. Since the polymer brush is stiff, it also imposes a 
minimum separation of 2h between particles, where h is the height of the brush. As a result, in the expression for 
eAB we can approximate the radial flexibility of the AB bond as 5 L — h. 
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FIG. 5: (Color Online) Comparison of the melting curves f{T) determined by our model to the experimental data of Chaikin 
et al(See Fig. 2 in [9]). The four data sets are for the four different polymer brushes used. For the model fits we find that 
(-^)2+~2.01 for crosses, 2.07 for solid triangles, 2.13 for empty triangles, and 2.35 for squares. 

We have now related the free energy eAB to the known thermodynamic parameters of DNA(AG = AH — TAS, 
AH = -77. 2-^ and AS = -227. 8;^^^), and the properties of linker DNA chains attached to the particles(grafting 
density cr ~ 3 x 10"^ linker length L ~ 20nm). The height of the polymer brush is /i = 13 ± 5nTO[9]. In 

fitting the experimental data we have taken the average value (h) = 13nm. Changing h within these bounds does 
not have a major effect on the melting curves. As a result there is one free parameter in the model, the confinement 
angle a. This angle determines (N) and AG, which in turn determine e^^s > and finally the melting profile /. 

With some minor modifications we can also analyze the "tail to tail" hybridization mode in a recent experiment 
of Mirkin et al[12]. In this experiment, R = 6.5nm gold nanoparticles were chemically functionalized with ss DNA 
linkers. The last 15 bases on the markers for particles of type A and B were chosen to be complementary to a 30 
base ss DNA linker. Since the strands are not ligated after hybridization, the experimental pictures are similar. 
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Temperature(°C) 



FIG. 6: (Color Online) The effect of the linker DNA grafting density a on the melting profile f{T). The results of the model 
are compared with experimental data in [12]. The three data sets represent grafting densities of 100% (squares), 50%(circles), 
and 33% (triangles) for which {N) ^,=2.Z2, 2.16, and 2.05 respectively. 



The unhybridized portion of the ss DNA linker simply serves as a spacer, and the hybridized portions become ds 
DNA, which we can again treat as rigid rods(see section II. B). This experiment is done without the addition of a 
polymer brush, but the grafting density is two orders of magnitude larger than the experiment of Chaikin et al. As 
a result, there is still an entropic repulsion[7] associated with compressing the particles below separation 2h. Here 
h could loosely be interpreted as the radius of gyration of the unhybridized portion of the linker. Despite the fact 
that L ~ i?, our planar calculation of AG provides a good fit to the experimental data. The other major difference 
is that now the attraction between particles is mediated by an additional DNA linker. 



AG = AGA + AGs-Tlog 



(57) 



The first term is the contribution to the free energy from the hybridization of the end sequence on linker A to the 
complementary portion of the 30 base ss linker. The hybridization free energies AGyi and AGs were calculated with 
the DINAMelt web server [20]. The last term is the contribution to the free energy from the translational entropy of 
the additional linker DNA, with cunk the additional linker concentration. This highlights some incorrect assumptions 
of the thermodynamic melting model[12], where the two hybridization free energies were not calculated separately, 
and the translational entropy of the additional linker DNA was ignored. By introducing dilutent strands to the 
system, one can probe the effect of the linker grafting density a on the melting properties of the assembly(See Figure 
2B in [12]). The agreement between the experimental data and our theory is good, except at small / values. This 
is not surprising, since comparing the two requires relating the measurement of optical extinction to the unbound 
fraction /. This is a nontrivial matter when dealing with aggregation, which corresponds to the small / regime. 



IV. CONCLUSION 



We have developed a statistical mechanical description of aggregation and melting in DNA-mediated colloidal 
systems. First we obtained a general result for two-particle binding energy in terms of DNA hybridization free 
energy AG, and two model-dependent parameters: the average number of available bridges (N) and the effective 
DNA concentration Ce//. We have also shown how these parameters can be calculated for a particular bridging 
scheme. In our discussion we have explicitly taken into account the partial ergodicity of the problem related to slow 
binding- unbinding dynamics. 
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In the second part of the paper it was demonstrated that the fractions of dimers, triniers and other clusters, 
inchiding the infinite aggregate, are universal functions of a parameter Zab/T = cab/T — log [ctot/2co]. We have 
applied the results of our theory to a particular scheme when the DNA bridge is made of two freely jointed rods (ds 
DNA). The obtained melting curves are in excellent agreement with two types of experiments, done with particles 
of nanometer and micron sizes. Furthermore, our analysis of the experimental data give an additional insight into 
microscopic physics of DNA bridging in these systems: it was shown that the experiments cannot be explained 
without introduction of angular localization of linker ds DNA. The corresponding localization angle a is the only 
fitting parameter of the model, which allows one to fit both the position and width of the observed melting curves. 

There arc several manifestations of the greater predictive power of our statistical mechanics approach, compared to 
the earlier more phenomenological models. First, once a is determined for a particular system, our theory allows one 
to calculate the melting behavior for an alternative choice of DNA linker sequences. Second, if the resulting clusters 
are separated, for example in a density gradient tube, the relative abundance of dimers, trimers, and tetramers can 
be compared to the values determined from the theory. 

Finally, the theory predicts aging of the colloidal structures, one experimental signature for which is hysteresis 
of the melting curves. Such an experiment proceeds by preparing a system above the melting temperature, and 
measuring the unbound fraction of colloids as the temperature is lowered. The system is allowed to remain in this 
cooled state for a very long time, perhaps months, during which multiple DNA bridges break and reform. During 
this time the colloids relax into a more favorable orientation state, including states which are not accessible by simply 
rotating about the contact point formed by the first DNA bridge between particles. This favorable orientation state 
is characterized by an average number of DNA bridges {N) greater than what we calculate in the partially crgodic 
regime. If the unbound fraction is then measured as the temperature is increased, the melting curve will shift to a 
higher temperature, consistent with a larger value of {N). 
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